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ABSTRACT 

We answer the long-standing question of which production mechanism is responsible for the cataclysmic 
variables (CVs) in globular clusters. Arguments have been given that range from mostly primordial presence 
to a significant contribution of later dynamical formation in close stellar encounters. We conclude, based on a 
thorough analysis of a homogeneous Chandra data set, that the majority of CVs in globulars has a dynamical 
origin. 

Subject headings: binaries: close — globular clusters: general — X-rays: binaries — cataclysmic variables 



1. INTRODUCTION 

In the 1970s, it was discovered that the number of outburst- 
ing low-mass X-ray binaries (LMXBs) per unit mass was or- 
ders of magnitude higher in globular cluster s than the rest 
of the Milky Way JClarkl Il975t lKat3 119751) . and theoreti- 
cally good reasons were quickly put forward to explain how 
LMXBs c ould be formed efficiently through encounters in 
globulars dFabian et al.ll975HHeggiell975HHiUsll975lll976t 
lSiitantvol1975r Taking an LMXB and replacing the neutron 
star by a white dwarf, we have a cataclysmic variable (CV). 
Whether globulars would also enhance CV formation through 
encounters was less clear. (Hut & Verbunt (1983) concluded 
that there still would be a c ontribution, but a much smaller 
one, and a detailed study by iDi Stefano & Rappaportl (119941) 
came to the same conclusion. 

There have been two main difficulties in testing this predic- 
tion. Until recently, globular cluster CVs have proven rather 
elusive, and only a handful were known in the entire Galac- 
tic globular cluster population. However, with high resolu- 
tion data from the Chandra X-ray Observatory and the Hub- 
ble Space T elescope, we are now finding dozens of CV s per 
cluster (e.g. lGrindlav et alll2001aHPoolev et aUl2002al) . An- 
other difficulty has been that the density of CVs in the field is 
not well known, so there is n o basis for comparison. Recently, 
iTownslev & Bildstenl (120051) have come up with a theoretical 
estimate for the number of CVs per unit mass of an old stel- 
lar population, and they concluded that the number of X-ray 
sources in 47 Tuc is compatible with their estimate, showing 
no overabundance relative to the field. This would mean that 
encounters have little to no effect on CV production, in con- 
trast to the outbursting LMXBs. 

On th e basis of Chandra observations of a dozen globular 
clusters JPoolev et al.H2003l paper I) found that the number of 
X-ray sources in a cluster scales with the encounter frequency 
(T) of the cluster significantly better than it scales with the 
mass of the cluster. This was the first clear proof that the 
population of X-ray sources in a globular cluster is largely 
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dynamically formed. 

Many of the X-ray sources in a number of clusters have 
been identified, and it has become clear that this is a hetero- 
geneous population, including LMXBs (in quiescence for the 
clusters under consideration here), CVs, chromospherically 
active main-sequence binaries, and millisecond pulsars. It 
was recognized that the quiescent LMXBs (qLMXBs) could 
be isolated from the rest of the population for the most part, 
and a number of authors noted the nea r-linear relationship 
between the number of qLMXBs and V JPoolev et al.ll2003t 
iGendre et al.ll2003l llMnke" et alJl2003aL . Because this rela- 
tionship was based on the number of total LMXBs (not just 
the outbursting ones), this firmly established that a higher rate 
of encounters resulted in a higher number of LMXBs, not just 
a higher fraction of outbursting systems, as could have been 
argued for the long-standing result from the 1970s. 

Having separated the qLMXB population from the general 
population, it is easily seen that qLMX Bs are the minority 
subpopulation. Because the correlation o flPoolev et al.lJ2003l) 
was based on the entire population, the implication is that, 
presumably, the other subpopulations must also receive some 
contribution to their formation from dynamical encounters. In 
this paper, we address the CV subpopulation. We sketch our 
approach in |2]and |3j analyse it critically in |4] and sum up 
in 33 

2. A CLEAR ENCOUNTER SIGNAL 

The main difficulty in our previous paper was the strong 
correlation between a cluster's mass (M) and its encounter 
frequency T, defined as J * {p 2 /v)Airr 2 dr, where r/, is the half- 
mass radius, p is the stellar density as a function of radius, and 
v is the velocity dispersion as a function of radius. This inte- 
gral is often approximated as po 2 r c 3 /vq, where r c is the core 
radius and the subscript zero refers to central values; this ex- 
pression simplifies to pa l5 r c 2 for a virialized system. By and 
large, clusters with a larger mass have a larger integrated col- 
lision number, and vice versa. Performing a correlati on test 
betwee n M and T for the 140 globular clusters in the iHarrisl 
( f!996l) catalog, we find a Spearman coefficient of 0.72 with a 
chance correlation probability of ~10~ 23 . However, there is a 
very large scatter in the correlation. At a given M, there are 
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TABLE 1 
Results of fitting n x = C x 



X 


c x 


X 2 /dof 


Q value 




I 


1.29 


27.8/5 


3.91 x 10" 


-5 


II 


3.99 


85.5/22 


1.97 x 10" 


-9 


III 


16.8 


69.1/18 


6.41 x 10" 


-8 


i+n+m 


22.0 


120.7/18 


1.12 x 10" 


-17 



orders of magnitude variation in Y; the best fit line between 
logM and logT has a mean absolute deviation of 0.69. 

In this paper, we alleviate the problem of the Y-M corre- 
lation by working in specific units. We define 7 = Y/M as 
the specific encounter frequency and n x = N x /M as the spe- 
cific number of sources of population x. The specific en- 
counter frequency 7 is a measure of the chance that a par- 
ticular star in a globular cluster undergoes an encounter, in 
contrast to Y, which is the chance that any star in a clus- 
ter has an encounter. We estimate the mass of the cluster as 
M/Mq = 3 x lQ0.MMv,@-Mv) us i n g t h e cluster a bsolute magni- 
tudes from the 2003 version of the catalog of lHarrisI dl996l) . 
and we have used units of 10 6 M Q to form our specific units. 

With these units, we can perform a very simple test. If pop- 
ulation x is primordial, then n x should be roughly the same 
for all clusters (where we neglect variations in the initial mass 
function (IMF), etc.; see 5j4J. This is easily checked by fit- 
ting for a constant, n x = C x . Should this give a statistically 
unacceptable fit, we conclude that population x is not simply 
primordial. To test for the effects of encounters, we can take 
the next step in complexity and assume some primordial con- 
tribution and some dynamical contribution: n x = C x +A/y ax . 

3. A (RELATIVELY) CLEAN SAMPLE OF CVS 

Determining how many of the Chandra X-ray sources in 
a globular cluster are CVs is difficult because most of the 
sources are discovered with only tens of X-ray counts. It is 
only through a comparison with HST data to look for blue 
or Ha bright sources within the X-ray error circles that most 
CVs are identified (this requires the matching of the Chandra 
and HST frames to sub-arcsecond precision, which is non- 
trivial). This is a time-consuming process and has only been 
accomplished for a large number of sources in just a few clus- 
ters. Here, we outline a method to obtain a relatively clean 
sample of CVs in each cluster on the basis of the X-ray data 
alone. 

In Fig. [0 we show an X-ray color-magnitude diagram for 
the X-ray sources within the half-mass radii of 23 globular 
clusters 1 . In total, ^600 sources are detected in these clus- 
ters, about 500 of which are represented in the figure (the rest 
are at lower L x ); we expect about 130 of these to be unrelated 
background sources b ased on the log /V- log S relationship of 
iGiacconi et alJ ( 1200 ll) . We denote those sources whose na- 
ture has been identified 2 by various symbols. We also plot 

1 NGCs 104 (47 Tuc), 288, 362, 5139 (w Cen), 5272 (M3), 5904 (M5), 
6093 (M80), 6121 (M4), 6205 (M13), 6218 (M12), 6266 (M62), 6341 (M92), 
6366, 6397, 6440, 6541, 6544, 6626 (M28), 6656 (M22), 6752, 6809 (M55), 

7099 (M30), and Terzan 5 

- The identifications come from 47 Tuc IGrindlav et all I2001at 
I Heinke et alJ [20051. NGC 288 IKong et all 12006ft. NGC 6093 (M80; 
IHeinke et all I2003hft. NGC 6121 ( M4: IBassa et all 2004). NGC 6397 
IGrindlav et all 12001 U). NGC 6440 lin't Zand et afl 120011) . NGC 6626 
M28: IBecker et alll2003ft . NGC 6752 IPoolev et aUI2002air and Terzan 5 
Wiina ndset all2005D . 



TABLE 2 
Results of fitting n x = C x +A x -y a ' 



X 


Cx 


a x 


X 2 /dof 


Q value 


I 
II 
III 
i+n+in 


40 +0 - 53 
27 +14 ° 

-12+ 14 

1 -235 

-5 3+" 


1 75 +0 - 43 
0.83«g 
0.27*g 
4S+ 017 


0.45/3 
27.7/20 
18.5/16 
22.4/16 


0.93 
0.12 
0.30 
0.13 



the locations of a number of field neutron star LMXBs in qui- 
escence. Further, we show the tracks of some spectral mod- 
els that have been successfully used to describe the spectra of 
quiescent neutron star LMXBs; these models consist of a ther- 
mal, blackbody-like component from the neutro n star surface 
propogated through a Hydgrogen atmosphere JZavlin et alJ 
1996) plus some contribution (up to about 50%) from a pow- 
erlaw component, whose origin is unknown. 

Based on observational considerations and the theoretical 
spectral models, we have divided the X-ray sources in glob- 
ular clusters into three populations, creatively denoted I, II, 
and III. Population I includes the two secure globular cluster 
qLMXBs (and most of the other known quiescent LMXBs 
from the field) and relatively little else. Population II in- 
cludes the rest of the sources down to a luminosity limit 
of 4 x 10 31 erg s" 1 (22 globulars have been observed to this 
limit), and population III goes a factor of ten lower to 4 x 
10 30 erg s" 1 (18 globulars have been observed to this limit). 
In population II, there have been a number of source iden- 
tifications. Although the constituents are mixed, the strong 
majority component appears to be CVs. This makes sense 
observationally, since few qLMXBs are as spectrally hard as 
most CVs, and the chromospherically active main-sequence 
binaries rare ly reach these luminosity levels dPempsev et alJ 
[T99llfl99l . There is little we can say about population III 
at this point except that it likely includes contributions from 
all known subpopulations. Differen ces in the X-ray lum inos- 
ity functions from cluster to cluster ( IPoolev et alJ2002bl) may 
suggest that the makeup of this population varies from cluster 
to cluster. 

4. ANALYSIS 

For each population x £ {I, II, III, I + 11 + 111} we plot n x as 
a function of Y (see Fig. |2j. We first fit a constant for each 
plot, which we find statistically unacceptable in all four cases. 
The goodness of fit and the likelihood of each fit accurately 
describing the data (the Q value, which is the probability of 
obtaining the observed x 2 assuming the model is correct and 
the best-fit parameters are the true parameters) are given in 
Table This can be interpreted as the likelihood of each of 
these populations being purely primordial. 

As discussed above, we can then fit n x = C V +A A 7 Q> , with the 
two terms representing contributions from a primordial pop- 
ulation (first term) and a dynamically influenced population 
(second term). In Table|2]we give the best fit parameters, the 
goodness of fit, and the probability that the model accurately 
describes the data. We plot the best fit relationship in Fig.|2] 

The re sults for »T+n+n are an extension of the work pre- 
sented in IPoolev et al.l ( 120031) . including more clusters, and 
working now in specific units. The results for the LMXB- 
dominated p opulation (wj) are al so an extension of the work 
presented in IPoolev et alJ ( 120031) and complementary to the 
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FIG. 1. — The Chandra color-magnitude diagram of 23 globular clusters, corrected for absorption. Above 4 X 10 30 erg s (4 X 10 31 erg s -1 ), there are ~500 
(~200) sources in this diagram, about 100 (15) of which are expected to be background sources. Not all clusters have been observed to the same limiting 
luminosity. Securely identified cluster sources are denoted by ■ for qLMXBs, A for CVs, for active main-sequence binaries, T for millisecond pulsars, and 
▼ for the pulsar in M28. In addition, we show the absorption-corrected locations of numerous field qLMXBs: • for pulsars and O for bursters. Tracks for NS 
atmosphere spectral models (plus some power law contributions) are also plotted; these spectral models have been used to successfully fit many qLMXBs. 



work o f lHeinkeetal1(l2006l) who also consider the effects of 
metallicity as a determining factor in the number of LMXBs 
present in a cluster; they find no evidence for metallicity de- 
pendence. 

The significant new result of this work is the demonstra- 
tion that the CV-dominated population («n) must also have 
a large dynamically-formed component. A primordial-only 
population is clearly ruled out, and our simple prescription for 
primordial-plus-dynamical is found to be acceptable. While 
our analysis has not taken into account a possibly confusing 
variation in IMF and other local primordial properties among 
globulars, such variations would not be expected to show the 
type of clear correlation between CV frequency and 7. 

A fair question to ask is whether population II is really dom- 
inated by CVs. This will only be resolved beyond a doubt 
when a significant majority of the sources in the population 
are identified. The other likely major constituent is active 
main-sequence binaries. As Fig.[T]indicates, few of the iden- 
tified main-sequence binaries in globular clusters are in this 
region, and, as noted above, very few active main-sequence 
binaries in the field reach luminosities above 4 x 10 31 erg s -1 . 
We therefore feel justified in identifying population II as CV- 
dominated. We do recognize that many of the non-CV sources 
in population II are binaries and, hence, can be expected to 
be overabundant because they also can be produced directly 
through encounters. 

Note that our conclusion about the dynamical formation of 
CVs is totally independent of the uncertain abundance of CVs 



in the field. It is therefor e neces sary to discuss the conclusion 
of ITownslev & Bildsten (2005J) that 47 Tuc shows no over- 
abundance of CVs. They present a prediction for the entire 
CV population, and here we deal with just the most X-ray 
luminous globular cluster CVs. Assuming a roughly simi- 
lar CV X-ray luminosity function in each cluster, this should 
not make a significant differ ence. Note that the prediction of 
ITownslev & Bildsten! J2005I) has a spread of a factor of three, 
and we are dealing here with numbers not much different than 
a factor of three. However, we have overwhelming statistical 
evidence that the trend we find is real, even though the results 
based on any individual cluster alone would not be meaning- 
ful at all. 

It may be th e cas e that the prediction of 
ITownslev & Bildsten! J2005I) just happened to be in the 
right range of 47 Tuc. For example, applying the prediction 
to to Cen, which is more than twice as massive than 47 Tuc, 
would imply twice as many CVs. In our CV-dominated 
population II, 47 Tuc has ~ 16 sources whereas 10 Cen has 
only ~2. The comparison of ITownslev & Bildsten! (120051) 
with 47 Tuc, while plausible for a single cluster, does not 
carry the statistical weight that we provide with our whole 
ensemble. 

5. SUMMARY 

We conclude that in dense globular clusters, such as NGC 
6388, NGC 6266 and 47 Tuc, the majority of CVs must be 
produced dynamically, by encounters between single stars 
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FIG. 2. — The specific number of sources for each population (%) plotted as a function of the specific encounter frequency (7), along with the best fit model 
of the form C x +A x -y a ' . a) Population I is almost exclusively quiescent LMXBs. The data have been binned for the model fits and for plotting, b) Population II 
is dominated by CVs. The data were binned for plotting, but the fits were performed on the unbinned data, c) Population III comprises many different source 
types. The red arrows indicate clusters which were n ot observed deeply enough to complete the census of this population, d) All sources above 4 X 10 30 erg s 
are shown, similar to Fig. 2 from Poolev et al. 1 2003) but here in specific units (i.e., dividing all values by cluster mass). 



and/or binaries. The obvious followup question is: what is 
the mechanism responsible? That will need to be sorted out 
in detailed simulations, in which each cluster is modeled sep- 
arately. This will be done, but it will take several years. How- 
ever, for now, we can already draw the firm conclusion that 
these CVs are predominantly formed dynamically by ruling 
out a primordial-only scenario, which has only a 1.97 x 10~ 9 
probability of describing the observations (TableQ. 

Because of the long-standing theoretical predictions that 
globulars should be rich in white dwarf binaries, they have 



been called Type la supernova factories iShara & Hurlevl 
120021) . We have now established observational evidence to 
support this idea. Globular clusters are adept at forming ex- 
otic binaries and their offspring: neutron star binaries, mil- 
lisecond pulsars, and, as we have now shown, white dwarf 
binaries. For this reason, globular clusters may play an in- 
teresting and sign ificant role in the story of th e short/hard 
gamma-ray bursts ( Zheng & Ramirez-Ruizl2006l) and provide 
exciti ng targets for gravitational wave detectors jBenacauistal 
[T999I 



REFERENCES 



Bassa, C, et al. 2004, ApJ, 609, 755 
Becker, W„ et al. 2003, ApJ, 594, 798 
Benacquista, M. 1999, ApJ, 520, 233 
Clark, G. W. 1975, ApJ, 199, L143 

Dempsey, R. C, Linsky, J. L., Fleming, T. A., & Schmitt, J. H. M. M. 1993, 
ApJS, 86, 599 

Dempsey, R. C, Linsky, J. L., Fleming, T. A., & Schmitt, J. H. M. M. 1997, 

ApJ, 478, 358 
Di Stefano, R., & Rappaport, S. 1994, ApJ, 423, 274 
Fabian, A. C, Pringle, J. E., & Rees, M. J. 1975, MNRAS, 172, 15P 
Gendre, B., Barret, D., & Webb, N. 2003, A&A, 403, LI 1 
Giacconi, R., et al. 2001, ApJ, 551, 624 

Grindlay, J. E., Heinke, C, Edmonds, P. D., & Murray, S. S. 2001, Science, 
292, 2290 



Grindlay, J. E., Heinke, C. O., Edmonds, P. D., Murray, S. S., & Cool, A. M. 

2001, ApJ, 563, L53 
Harris, W. E. 1996, AJ, 112, 1487 
Heggie, D. C. 1975, MNRAS, 173, 729 

Heinke, C. O., Grindlay, J. E., Lugger, P. M., Cohn, H. N„ Edmonds, P. D., 

Lloyd, D. A., & Cool, A. M. 2003, ApJ, 598, 501 
Heinke, C. O., Grindlay, J. E., Edmonds, P. D., Lloyd, D. A., Murray, S. S., 

Cohn, H. N„ & Lugger, P. M. 2003, ApJ, 598, 516 
Heinke, C. O., Grindlay, J. E., Edmonds, P. D., Cohn, H. N„ Lugger, P. M., 

Camilo, F, Bogdanov, S., & Freire, P. C. 2005, ApJ, 625, 796 
Heinke, C. O., Wijnands, R., Cohn, H. N., Lugger, P. M., Grindlay, J. E., 

Pooley, D., & Lewin, W. H. G. 2003, ApJ, submitted 
Hills, J. G. 1975, AJ, 80, 809 
Hills, J. G. 1976, MNRAS, 175, IP 
Hut, P., & Verbunt, F. 1983, Nature, 301, 587 



Dynamical Formation of Globular Cluster CVs 



5 



in't Zand, J. J. M., van Kerkwijk, M. H., Pooley, D., Verbunt, F., Wijnands, 

R., & Lewin, W. H. G. 2001, ApJ, 563, L41 
Katz, J. I. 1975, Nature, 253, 698 

Kong, A. K. H., Bassa, C, Pooley, D., Lewin, W. H. G., Homer, L., Verbunt, 
R, Anderson, S. F., & Margon, B. 2006, ApJ, submitted, astro-ph/0603374 
Pooley, D., et al. 2002a, ApJ, 569, 405 
Pooley, D., et al. 2002b, ApJ, 573, 184 
Pooley, D., et al. 2003, ApJ, 591, L131 (paper I) 
Shara, M. M., & Hurley, J. R. 2002, ApJ, 571, 830 



Sutantyo, W. 1975, A&A, 44, 227 

Townsley, D. M., & Bildsten, L. 2005, ApJ, 628, 395 

Wijnands, R., Heinke, C O., Pooley, D., Edmonds, P. D., Lewin, W. H. G, 

Grindlay, J. E., Jonker, P. G, & Miller, J. M. 2005, ApJ, 618, 883 
Zavlin, V. E., Pavlov, G. G, & Shibanov, Y. A. 1996, A &A, 315, 141 
Zheng, Z., & Ramirez-Ruiz, E. 2006, ApJ, submitted, astro-ph/0601622 



